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Abstract 

We study magnetic field effects on spontaneous Fermi surface symmetry breaking with d-wave 
symmetry, tlie so-called d-wave "Pomeranchuk instability" . We use a mean-field model of electrons 
with a pure forward scattering interaction on a square lattice. When either the majority or the 
minority spin band is tuned close to the van Hove filling by a magnetic field, the Fermi surface 
symmetry breaking occurs in both bands, but with a different magnitude of the order parameter. 
The transition is typically of second order at high temperature and changes to first order at low 
temperature; the end points of the second order line are tricritical points. This qualitative picture 
does not change even in the limit of a large magnetic field, although the magnetic field substantially 
suppresses the transition temperature at the van Hove filling. The field produces neither a quantum 
critical point nor a quantum critical end point in our model. In the weak coupling limit, typical 
quantities characterizing the phase diagram have a field-independent single energy scale while its 
dimensionless coefficient varies with the field. The field-induced Fermi surface symmetry breaking 
is a promising scenario for the bilayer ruthenate Sr3Ru207, and future issues are discussed to 
establish such a scenario. 

PACS numbers: 71.18.+y, 75.30.Kz, 74.70.Pq, 71.10.Fd 
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I. INTRODUCTION 



Usually the Fermi surface (FS) respects the point-group symmetry of the underlying lat- 
tice structure. However, recently a symmetry-breaking Fermi surface deformation with a d- 
wave order parameter, where the FS expands along the direction and shrinks along the ky 
direction, or vice versa, was discussed in various two-dimensional interacting electron model 



on a square lattice, the t-J 



|l, 2, 3] Hubbard, Q,H, 3 



and extended Hubbard [7] model. This 



d-wave type Fermi surface deformation ((iFSD) is often called (i-wave Pomeranchuk instabil 



ity, referring to Pomeranchuk's stability criterion for isotropic Fermi liquids. 



Bl] However, the 



Fermi surface symmetry breaking can happen even without breaking such a criterion, since 
the instability is usually of first order at low temperature, jo], lo| Moreover, the new concept 
of the Fermi surface symmetry breaking is applicable also to strongly correlated electron 
systems such as those described by the t-J mo The (iFSD instability is driven 

by forward scattering processes of electrons close to van Hove points in the two-dimensional 
Brillouin zone. The instability is thus purely electronic and the lattice does not play a role. 
As a result, symmetry of an electronic state is reduced from 6*4^ to C2V, while the lattice 
retains symmetry as long as no electron-phonon coupling is considered. 
The d 



and t- 
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''SD competes with superconductivity. Several analyses of the Hubbard 
3j model showed that the rf-wave superconductivity becomes a leading instability 



and the spont 



aneous c/FSD does not happen. However, appreciable correlations of the 



c/FSD remain. ISj As a result, the system becomes very sensitive to a small external xy- 



anisotropy and shows a giant response to it, leading to a strongly deformed FS. This idea 
was invoked for high-temperature cuprate superconductors. In particular, the recently 
observed anisotropy of magnetic excitations in YBaoCu.sOfi+.x 1 14] has been well understood 
in terms of c/FSD correlations in the t-J model. ISj Although the reduced symmetry due 
to the (iFSD is the same as the electronic nematic phase proposed in the context of the 
so-called spin-charge stripes, [l^ the underlying physics is very different. The Fermi surface 
symmetry breaking does not require the assumption of charge stripes, but is driven by 
forward scattering processes of electrons. The dFSD provides an essentially different route 
to the nematic phase. 

Although a spontaneous rfFSD has not been proposed for cuprates because of the com- 
petition with the (i-wave singlet pairing, the material Sr3Ru207 turned out to have an 
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interesting possibility of a spontaneous c?FSD.[17| Sr3Ru207 is the bilayered ruthenate with 
two metalhc Ru02 planes where Ru-ions form a square lattice. Unlike the corresponding 
single-layered material Sr2Ru04, a well-known triplet superconductor, jisl Sr3Ru207 has a 
paramagnetic ground state. 1^ However, the material is close to a ferromagnetic transi- 
tion, which was su gge sted by the strongly enhanced uniform magnetic susceptibility with a 
large Wilson ratio, 20l | uniaxial-pressure-induced ferromagnetic transition, 2l|] inelastic neu- 



tron scattering. 
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2J] By applying a magnetic field h, 



22l | and band structure calculations. 
Sr3Ru207 shows a metamagnetic transition aXh = around which non- Fermi liquid behav- 
ior was observed in various quantities: resistivity, [25|, l26| specific heat,|26|, |27|, l28| thermal 



expansion. 



29( 1 and nuclear spin-lattice relaxation rate. 30|| This non-Fermi liquid behavior 



was frequently discussed in terms of a putative metamagnetic quantum critical end point 
(QCEP), and in fact Sr3Ru207 was often referred to as a system with a metamagnetic 
QCEP.H However, after improving sample quality, the hypothetical QCEP turned out to 
be hidden by a dome-shaped transition line of some ordered phase around /ic-[17|] While 
a second order transition was speculated to occur around the center of the dome, a first 
order transition was confirmed at the edges of the transition line and was accompanied by a 
metamagnetic transition. Grigera et al. 17| associated this instability with the spontaneous 
c/FSD, which turned out to be consistent with a large magnetoresistive anisotropy recently 
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observed inside the dome-shaped transition line. 

Any first order transition as a function of a magnetic field is generically accompanied by a 
metamagnetic transition, which follows from the stability of the thermodynamic potential. 
This was demonstrated in the case of the first order dFSD transition in connection with 
Sr3Ru207.[33| Quite recently we showed that the most salient features observed in Sr3Ru207, 
not only the metamagnetic transition but also the phase diagram and the non-Fermi liquid 
like behavior of the uniform magnetic susceptibility and the specific heat coefficient, are well 
captured in terms of the dFSD instability near the van Hove singularity without invoking 
a putative QCEP. 34] We also predicted anomalies associated with the ciFSD instability in 
the temperature dependence of the magnetic susceptibility and the specific heat. 34] 

The main purpose of this paper is to ex pan d our previous paper about the rfFSD in- 
stability in the presence of a magnetic field. 3J] It is particularly interesting to perform a 
comprehensive analysis of magnetic field effects on the dFSD instability, since a magnetic 
field is often employed as a tuning parameter of a quantum phase transition. A naive ques- 
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tion may be whether a QCEP is reahzed for the liFSD instabihty as for the ferromagnetic 



instabihty. 
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361] The phase diagram of the c/FSD is known to be characterized by several 
universal numbers, lo| which can be compared directly with experimental data. It is then 
interesting also how the universal numbers evolve in the presence of a magnetic field. 

We analyze the (iFSD instability in the charge channel in a one-band model on a square 
lattice. The model describes electrons interacting via a pure forward scattering interaction 
driving the (iFSD in the presence of a magnetic field (Sec. 11). We solve this model numeri- 
cally in Sec. Ill and investigate the weak coupling limit analytically in Sec. IV. In Sec. V, we 
discuss the reportec phase diagram for Sr^RuoO^jl?] as well as relations to other scenarios 
such as a QCEP, 3l|, |35|, |37t| phase separation, [38] and magnetic domain formation. 39| Sec- 



tion VI is the conclusion. In Appendix A, detailed features of the ciFSD phase diagram are 
presented. In Appendix B, we analyze the (iFSD instability in the spin channel, often called 
spin- dependent Pomeranchuk instability, which shows exactly the same phase diagram as in 
the charge channel. 



II. MODEL AND FORMALISM 



We investigate the (iFSD instability in the charge channel under a magnetic field on a 
square lattice. The minimal model reads 

-f^ = - Ai)^ka + ^ ^ /kk' ^ka?^k'a' - ^ O-rika (l) 

k,cr k,cr,k',o-' k,(7 

where riko- = cj^.^Cko- counts the electron number with momentum k and spin a; cj^^ (ck<j) is 
an electron creation (annihilation) operator; fi is the chemical potential; is the number 
of lattice sites; h is an effective magnetic field and is defined as h = ^qhbH where g is a 
(7-factor, /i^ is Bohr magneton, and if is a magnetic field. For hopping amplitudes t and 
t' between nearest and next-nearest neighbors on the square lattice, respectively, the bare 
dispersion relation is given by 

ek = — 2t(cos kx + cos ky) — At' cos k^ cos ky . (2) 

The forward scattering interaction driving the spontaneous (iFSD has the form 

/kk' = -g dkd]^' , (3) 
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with a coupling constant g > and a d-wave form factor (ik = cos — cos ky. This ansatz 
mimics the structure of the effective interaction in the forward scattering channel as obtained 
fo. .Ke ...Q Hubb.d.Q and extended Hnbb.dQ .odel. Fo. . ^ 0. th. .ode, and a 
similar model were studied in Refs. HQ, and 9,, respectively. 

We decouple the interaction by introducing a spin-dependent mean field 

= -^J^^Ka), (4) 

k 

which becomes finite when the system breaks orientational symmetry and is thus the order 
parameter of the rfFSD. The mean-field Hamiltonian reads 

TV 

^^MF = ^ka n]^a + — , (5) 

k,cr ^ 

where 

^ka = ^k + Vdk - yu<7 • (6) 

Here the a-summed mean filed 1] = rj^j enters ^ko-, and thus a finite r/o- in general induces 
a finite ?7_o-; the magnetic field is absorbed completely in the effective chemical potential 
yUo- = /i + crh. The grand canonical potential per lattice site is given by 

'^ = -^EMl + e-^-/^) + |^. (7) 

kcr ^ 

By minimizing Eq. ([7]) with respect to rj, we obtain a self-consistency equation 

k,cr 

We consider the solution with rj > 0, since the free energy Eq. ([7]) is an even function with 
respect to rj. The self-consistency equation is written also as 

^'^ = -f (9) 

k 

Note that our Hamiltonian ([1]) does not allow momentum transfer, and thus the mean-field 
theory solves our model exactly in the thermodynamic limit. 



III. NUMERICAL RESULTS 



LDA band calculations 23|, |2J] for Sr3Ru207 without a magnetic field showed that the 
electronic structure is similar to that for the single-layered material Sr2Ru04 except that 
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there are six FSs because of the bilayered structure. Since the dFSD instabihty is driven 
by electrons near the van Hove points on a square lattice, we focus on the FS closest to 
k = (vr, 0) and (0,7r), and mimic such a FS by choosing t' /t = 0.35. For h = the bare 
dispersion has the van Hove energy at At' = lAt, from which we measure the chemical 
potential fi. We take g/t = 1 for numerical convenience, but the result for g/t = 0.5 shall 
be mentioned in the context of Fig. HI We set t = 1 so that all quantities with dimension 
of energy are in units of t, and consider a region of > in this paper since the result is 
symmetric with respect to h ^ —h and cr —a. 

Figure [D^a) shows a phase diagram for ^ = —0.4 in the plane of the magnetic field h 
and temperature T; the dotted line at h = h^n = 0.4 represents the van Hove energy of 
the up-spin band (majority band). The dFSD transition is of second order for high T and 
changes to first order for low T; the end points of the second order line are tricritical points. 
Hence the transition line forms a domed shape around the van Hove energy. Figure [U^b) 
shows the h dependence of the order parameter t] for low T together with r]^ and ?7|. We 
see that although the phase transition happens around the van Hove energy of the up-spin 
band, both rj^ and rji show a jump at the first order point; the ?7| has the same sign as rji, but 
with a different magnitude. The FSs at T = 0.01 are shown in Fig. [T]^c) and (d) for h = 0.3 
and 0.5, respectively, together with the FSs for g = 0; the outer (inner) FS corresponds to 
the up-spin (down-spin) electron band; the splitting of the FSs is due to the Zeeman energy. 
The FS instability drives a deformation of both FSs and typically leads to an open outer 
FS. Electron density for each spin also shows a jump at the first order transition point, but 
the size of the jump is different [Fig. [I^e)]. This difference yields a metamagnetic transition 
as shown in Fig. This is a generic consequence of the concavity of the grand canonical 
potential when a first order transition occurs as a function of h, the canonical conjugate 
quantity to the magnetization m. The metamagnetic transition disappears for high T where 
the (iFSD transition becomes second order. The T dependence of the magnetization is shown 
in Fig. [11(g) for several choices of h. The magnetization shows a kink at the second order 
transition. Roughly speaking, compared with the non-interacting case, the magnetization is 
enhanced (suppressed) by the ciFSD transition for h < h^u{h > Ku). Since the transition 
is of first order for low T as a function of h, there appears a phase separation as a function 
of the magnetization [Fig. [11(h)]. The width of the phase separated region corresponds to a 
magnitude of a jump of m at the first order transition point seen in Fig. [11(f). 
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The (iFSD instability has a dome-shaped transition hne around the van Hove energy of 
the up-spin band for the chemical potential {fi < 0) as shown in Fig. [TJ^a). When one invokes 
much larger yu(> 0), the dFSD transition then happens around the van Hove energy of the 
down-spin band. When fi is around zero, the second order line extends down to h = and 
a first order line appears only on the larger h side for low T as shown in Fig. [2J For some 
special values of /i, see Appendix A. 

Although a metamagnetic transition generically accompanies a first order (iFSD phase 
transition, the metamagnetic transition can also occur inside the symmetry-broken phase 
because of a level crossing between two local minima of the free energy. This is actually the 
case in the present model as demonstrated in Fig. [31 where m is plotted as a function of h 
for fi = —0.2; the h dependence of the order parameter is shown in the inset. In addition 
to a jump associated with the first order c/FSD transition a.t h ^ 0.34 (see also Fig. [2]), 
another metamagnetic transition appears at /i ~ 0.28 in the symmetry-broken phase. This 
metamagnetic transition is a weak feature in the sense that it is smeared out by thermal 
broadening for a relatively low temperature while the metamagnetic transition associated 
with the (iFSD instability from the symmetric state is robust up to T = T*". 

Figure Hl^a) summarizes numerical results in the plane of {fi, h) at low T = 0.01; the 
solid circles denote a first order transition to the symmetry-broken phase, where the order 
parameter t] shows a jump from zero to a finite value, accompanied by a metamagnetic 
transition as a function of h; the cross symbols are positions, where a first order transition 
occurs in the symmetry-broken phase due to a level crossing of local minima of the free 
energy, yielding an additional metamagnetic transition; the dotted line represents positions 
of the van Hove energy for the up-spin band (/i < 0) and the down-spin band (^u > 0). We 
see that the dFSD phase is stabilized around the van Hove energy of each spin band as a 
function of h for a given fi. With respect to the axis /i = 0, the phase diagram is nearly 
symmetric and becomes fully symmetric when t' is set to be zero. When we take a smaller g, 
the symmetry-broken phase is stabilized closer to the van Hove energy, but the qualitative 
features of Fig. Hl^a) do not change except that the energy scale is substantially decreased. 

We define Tj" as at the van Hove energy, namely a.t h = \fi\, and show its h dependence 
in Figs.m^b) and (c) for /i < and /i > 0, respectively. The TJ^ is suppressed with increasing 
magnetic field. In particular, for a smaller g, a relatively small h suppresses TJ^ drastically. 
However, the suppression saturates for larger h, leading to a finite TJ^, where we obtain a 



7 



phase diagram similar to Fig. [T]^a). That is, neither a quantum critical point (QCP) nor a 
QCEP of the (iFSD is realized by the magnetic field, which we will further discuss in Sec. IV. 



IV. WEAK COUPLING LIMIT 

The (iFSD instability occurs around the van Hove filling and thus the transition is dom- 
inated by states with momentum near the saddle points of e^. In the weak coupling limit, 
therefore, the mean-field equations can be treated analytically by focusing on the state near 
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in the absence of a magnetic field. Since 



the saddle points, similar to the analysis in Ref . 
the magnetic field is absorbed completely in the effective chemical potential fi^ = fJ^ + crh, we 
can extend such an analysis to the present case by allowing the a dependence of the chemical 
potential. The chemical potential fi is measured from the van Hove energy at /i = so that 
/io- = indicates that the a-spin band is at the van Hove filling. We first determine a zero 
temperature phase diagram in the plane of {^,h), and then investigate Tc suppression by 
a magnetic field, /i dependence of typical quantities characterizing the phase diagram, and 
the limit h ^ oo. 



A. Zero temperature phase diagram 



Following the analysis in Ref 



id, the self-consistency equation Eq. ([8]) is written as 



V 



? [it^^ - V) log IfJ'a -v\- if^a + v) log l/ia + r^l + 2?] (1 + log ] 



(10) 



where g = 2mg/Ti'^ is the dimensionless coupling and ca = A^/(2m) is a cutoff energy; m is 
the effective mass near the van Hove energy and is related to the hopping integrals t and t' . 
The grand canonical potential is then given by 

" 1 1 

2^" 2 

'1 



uj{r];fi,h) 



2m 



- (1 + logeA) 



4E 



2 (^'x + rjf log l/ia + + vf^og - 7]\ 



+ const. (11) 



It is not difficult to see that Eqs. (fTOl) and ( fTTl) are symmetric with respect to interchange of 
h and /i, that is, the magnetic field h plays exactly the same role as the chemical potential /x. 
We focus on the region < < h and introduce rescaled variables fj = rj/h and fi^ = jl + a 
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with fi = fi/h. Equation (fTOj) then reads 

^ = \Y1 [(^'^ ~ ^) lAa - ^1 - (/i<x + ^) log l/i^ + vl] (12) 
with a renormalized coupling constant 

r'=9-' + '2\ogh-2{l + log eA) . (13) 
Similarly Eq. ( ITTi) is written as Cij(r7; /i, /i) = ^h'^oj^f]; fi), where 



1 ~ ~ 2 ~ 1 2 ~ 

- {fia + 7]) log l/i^ + r/l + - (/i^ - r]) log l/i^ - f]\ 



+ 9E 



(14) 



and the energy is shifted such that u;(f/ = 0; /i) = 0. 

At zero temperature, the rfFSD transition is usually of first order as we have seen in 
Figs. [T](a) and[2l The first order transition is determined by solving Eq. (fT2|) and Cj{fi] /i) = 
numerically for a given /i, yielding a solution rfi and gi. A corresponding magnetic field hi 
is then obtained from Eq. fll3p 

hi = exp (l + — ^-^^ eAe-i/(^^) , (15) 



2^1 (/i) 

and thus the chemical potential and the order parameter are = jlhi and r/i = rfi/ii, 
respectively. All quantities, hi, /xi, and 771, are determined by a single energy scale cac"^/*^^^^ 
and the magnetic field just changes its dimensionless coefficient. We plot (/xi, hi) in Fig. 
the phase diagram is symmetric with respect to the /i = axis and the h = axis; 
such symmetry is not seen in our numerical result [Fig. Hl^a)], since a finite t' breaks the 
symmetry when g is not in the weak coupling limit. The first order transition line has a 
kink at ^ LScac"^/'^^^^ and h ^ 0.54eAe~^/'^^^^ (see the inset of Fig. [5]). At this point, the 
solution of Oj{fj; fi) = [Eq. ( IT^ ] shows a jump, indicating double local minima of the free 
energy oj{fj) away from fj = 0, which then may yield an additional first order transition in 
the symmetry-broken phase as seen in Figs. [3]and|l](a). 



B. Tc suppression by a magnetic field 



A magnetic field suppresses the (iFSD transition temperature [Fig. Hl^b)]. Here we clarify 
key factors of this suppression. 
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Since the dFSD transition is of second order as a function of T at the van Hove fiUing, 
the Tc is obtained by hnearizing the right-hand side of Eq. ([8]) with respect to i], namely 



l = (7iV2(/i,/i,Te), (16) 

where we introduce 

iV,(/i, /z, T) = -^Y. - /i.) , (17) 

k;,cr 

a weighted density of states averaged over an energy interval of order T around p is an 
even integer; /' is a first derivative of Fermi distribution function with respect to e^. In the 
weak coupling limit [10], Eq. (IT7|) reads 

No p dependence appears on the right-hand side, since we have redefined = ^{coskx — 
cos ky) in the present analysis in the weak coupling limit[10] so that \d]^\ = 1 at k = (vr,0) 
and (0,7r), namely at the saddle points of e^- The van Hove filling of the a spin electron 
band is set by choosing 

Ha = fJ^ + o'h = , = 11 — ah = —2ah . (19) 

We consider the a =| case, namely for /x > and h > 0. Since 

- rdelog^/'(e-2/i)^^^^^logeA + a(/i,T), (20) 



where 



we obtain 



roo d 1 

a{h,T) = J dx\og\Tx + 2h\ — -^, (21) 



^2 



2\og-^ + a{h,T)-a{0,T) 



(22) 



where a(0,T) = log[2e'^/(7rT)] with Euler constant 7 ~ 0.577. Defining a rescaled variable 
hx = h/T, we may write 



C(M = a{0,T) - a{h,T) 
2e^ _ r 



2e^ f°° ~ d 1 

log / dxlog|x + 2/iT|^— — , (23) 
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where Cihr) increases monotonically with and C(0) = and C(^r) = C(~^t)- Substi- 
tuting Eq. (l22l) into Eq. (fT6|l . we obtain TJ^ at the van Hove fiUing for a given hx 

t;h = '^,^^-vm)^-\ahT) (24) 

TT 

= T;«(0)e-5«''^) (25) 

where 

T;H(o) = ^e^Q-Vm (26) 

TT 

is the critical temperature for = 0; the corresponding magnetic field is obtained as /i = 
hxTJ^. We evaluate C(^t) numerically and show in Fig. [6] the h dependence of TJ^. The 
TJ^ is suppressed with as we have seen in Fig. 111(b). Defining a magnetic field /ii/2, at 
which TJ^ is suppressed down to a half value of TJ^(O), we read off /ii/2 = l.OOTJ^(O) from 
Fig. [6l Therefore from Eq. fl26|) we obtain 

oc eAe-^/(25) , (27) 

that is, the suppression of T^^ is controlled by the dimensionless coupling constant g = 
2mg/iT'^. When g becomes smaller, hi/2 gets smaller exponentially, leading to a strong 
suppression of TJ^ by the magnetic field. Similarly, the smaller effective mass suppresses Tc 
substantially with a magnetic field. 



C. /X dependence of characteristic quantities of the phase diagram 

We have derived the analytic expressions for hi [Eq. ([15])] and T^^ [Eq. i^]. There 
is another important quantity characterizing the dFSD phase diagram, the tricritical point 
{T^^\ htri), which we first calculate for a given fi. We then summarize how hi, TJ^, T^^\ and 
htvi evolve as a function of /i. Since all these quantities are scaled by a single energy, ratios 
of these different quantities become universal, whose fi dependence is also clarified. 

At the tricritical point, both the quadratic and quartic coefficient of the Landau free 
energy, uj{ri), vanish simultaneously. This condition is determined by Eq. f|T6l) and 

£^N,ifi,h,T) = (28) 

where A^4 is defined in Eq. f fTTj) . While we have analyzed Eq. f|T6|) at the van Hove point 
[Eq. (fT9l) ] in the previous section, we here consider Eq. (fT6l) for general /ig- and obtain 

T*" = e"eAe-^/(2^) ^ (29) 
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where 

a{P'T,hT) = J dx log \x + fiT + ^^t\^ ^^ ^ ^ (30) 
with fiT = yu/T". Similarly Eq. fl25]) is written in the weak coupling limit as 

It is easy to observe that Eqs. fl30l) and fl3T]) are symmetric with respect to fir ^ /^r and 
thus h plays exactly the same role as fi. We solve Eq. (13T|) numerically for a given /j-t and 
determine /i^- The tricritical temperature T*" is then obtained from Eqs. fl29|) and fl30l) : 
the original yU and h are determined as /itri = hxT^^^ for a given = JitTI^^. Figure Cl^a) 
summarizes TJ^, T*", /itri, and hi as a function of /i. We see that all quantities are suppressed 
with increasing but do not reach zero; in fact they approach certain asymptotic values 
as we show in the next subsection. Since all quantities are scaled by the single energy, we 
obtain universal numbers by taking ratios of the different quantities. In Fig. [7](b), we plot 
representative universal ratios Tl^^jT^^ and T*"/|/itri — /ivnl, where /Ivh is the van Hove 
energy and is given by /ivH = The strong /x dependence appears for relatively small ix 
and the variation is within a factor of 1.5. 

Although universal ratios are obtained in the weak coupling limit, these numbers also 
characterize well the phase diagram for a relatively large g. For example in Fig. [T](a), 
r*"/|/itri — ^vh| ~ 0.5 — 0.6 and T*"/T^^^ k, 0.5 — 0.75, which are comparable values to those 
in Fig.EKb). 



D. Large h limit 

For a larger /i, each quantity in Fig. [7](a) approaches a certain asymptotic value. In 
addition, the universal ratios in Fig. Cl^b) converge to the same values as those at /i = 0. To 
understand such asymptotic behavior, we consider TJ^ in the limit /i| ^ oo while keeping 
/i| = [see Eq. f|T9|) ]. namely, the limit oi ^ oo and h ^ oo. Since the up-spin electron 
band becomes fully occupied, we have /'(e — /i^) = for |e| < ca. Hence from Eq. (ITSl) we 
have 

N2{T) f\e\ogf-f{e) (32) 

^ J —eh 1^1 



eA/T^oo 2m 2eAe'^ 
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The gap equation Eq. (fT6ll then yields 

= — eAe-i/5 (34) 

TT 

= r;«(0)e-i/(2§) _ (35) 

That is, the (iFSD transition occurs even for h ^ oo under the condition of /i = |/i|, although 
TJ^ is suppressed a factor of e"^/^^^-* compared to the case of h = 0. Since the other spin 
band is fully occupied (empty) for fi +oo (/x —oo), only one spin band is subject to the 
ciFSD instability. Therefore in the h ^ oo limit, our model is reduced to a "spinless" fermion 



model, and thus the same results as those foTh = (Sec. V in Ref. lid) are obtained except 
that the energy scale is replaced by e^e'^^^. Hence various ratios of different quantities 
characterizing the phase diagram show exactly the same universal numbers as those for 
h = 0. The magnetic field just reduces the energy scale and cannot produce a QCP nor a 
QCEP of the dFSD instability. 



V. RELEVANCE TO SraRuaOr 

We have shown that when either the up- or the down-spin electron band is tuned by a 
magnetic field close to the van Hove filling, the rfFSD transition occurs in both bands but 
with a different magnitude of the order parameter. The field-induced dFSD instability is a 
promising scenario for Sr3Ru207. In fact, several important features observed experimentally 
such as the phase diagram, the metamagnetic transition, and the non-Fermi liquid like 
behavior of the magnetic susceptibility and the specific heat coefficient are well captured in 
terms of the rfFSD instability around the van Hove filling as we have shown in Ref. 3^. 

On the basis of the present results, we discuss in more detail the reported phase diagram 

n 

for Sr3Ru207.p/7| Since no experimental evidence of a symmetry-broken phase was obtained 



at h = and LDA band calculations 
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2J] showed that the van Hove energy is located 



above the Fermi energy, corresponding to /i < 0, we expect that the chemical potential n 
is away from the van Hove energy, but rather close to it, for example fi ~ —0.4 for the 
parameters shown in Fig. |H Figure [T](a) is a representative phase diagram of the rnagnetic 



field-induced rfFSD instability, which is very similar to the reported phase diagram. 17| The 
maximal Tc in the experiment is about 1 K, which is much smaller than the energy scale in 
Fig.[l](a). The coupling constant g in Sr3Ru207 is thus expected to be very small. In the weak 
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coupling limit (Sec. IV), the phase diagram is characterized by a single energy scale eAe~^/^^ 
and thus various ratios among TJ^, T*''', htri, and hi become universal numbers [Fig. Cl^b)]. 
The universal ratios are compared directly with experimental data. The data by Grigera et 
a/.Q provide T^^ ^ 1 K, T*" ~ 0.6 K, ht^ ^ hi = \q^ibH with H ^ 7.8 and 8.1 Tesla, 
and /ivH = Ifl/iB^^vH with H^n ~ 7.95 Tesla. We thus obtain ht.Jhi ^ 1, Tl''/T^^ ~ 0.6, 
and T*"/|/itri — h^n\ ~ (^b • 0.6)/(|g;UB ■ 0.15) ~ 12g~^ with Boltzmann constant ks- The 
first one is consistent with our result Fig. Cl^a); the value of the second one is comparable 
with Fig. [71(b); as for the last one, however, the discrepancy is by a factor of 10 if we assume 
= 2. 

While the field-induced c/FSD instability is well analyzed in the present model, further 
studies are necessary to make more quantitative comparison with experiments, (i) The 
interaction in our model retains SU(2) symmetry and thus the present theory cannot address 
the issue why the possibility of the dFSD instability was clearly observed only over a narrow 
region of applied field angle to the Ru02 plane. j4o| One possible origin of such a field 
angle dependence lies in a magnetic anisotropy, which we expect to originate mainly from a 
relatively strong spin-orbit coupling connected with the heavy Ru-ion. Our model should be 
extended by taking the magnetic anisotropy into account for a more quantitative study, (ii) 
Inclusion of a magnetic interaction in our model rnay also be necessary, since Sr3Ru207 is 



mm 



24j This is suggested also 



expected to be close to a ferromagnetic transition. |20, 

by comparing Fig. [11(g) with the experimental data. 17| In Fig. [11(g) the dFSD instability 
produces a kink in the T dependence of m at the transition temperature T^. This kink 
is actually observed in the experiment IT]] and our result for h > 0.35 is similar to the 
experiment. But the present theorycannot reproduce the observed upward curvature of 
the T dependence of m for T > Tc, [17] which may come from a magnetic interaction, (iii) 
Sr3Ru207 is a bilayered material. Since the (iFSD instability is driven by forward scattering 
processes of quasi-particles near k = (tt, 0) and (0,7r), a form of bilayer coupling can be 
important if the bilayer coupling gives rise to dispersion around k = (tt, 0) and (0, tt). 



Insights into the kz dispersion will be obtained from further detailed LDA calculations. 



J 



le dFSD instabi 
a y Hubbard 



ity was discussed in basic lattice models such as two-dimensional t- 



[ty was aiscussea m Dasic latti 
y, and extended Hubbard [3] 



and extended Hubbard [7|] model, and can be a generic tendency 
near van Hove filling in correlated electron systems. j4l| Hence the c/FSD is an interesting 
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possibility when the Fermi energy is tuned close to the van Hove energy in other materials 
also such as Sr2_yLayRu04 42| where La-substitution introduces electron carriers and makes 
the FS closer to the van Hove point. However, La introduces some disorder in the RUO2 
plane; its effects should be considered carefully, since the physics near the van Hove singu- 
larity may in general depend strongly on sample purity. In fact, the specific heat coefficient 



for low T in Sr2_yLayRu04 



42| shows different behavior from that in Sr3Ru207 although 



both systems are expected to be nicely tuned close to the van Hove filling. It is interest- 
ing to investigate impurity effects on the (iFSD instability by chemical (electron) doping to 
Sr3Ru207. According to our phase diagram [Figs. Il](a) and [5], see also Appendix A], the 
c/FSD instability can occur for a smaller magnetic field and even without the field if the 
impurity effect by the electron doping is not serious. 

Around the van Hove filling, various ordering tendencies such as antiferromagnetism. 



ferromagnetism, superconductivity, and (i-density wave develop. [111. Il2l| and compete with 
the (iFSD instability. Since the (iFSD instability is suppressed by a magnetic field and its 
suppression is controlled by ^ oc mg [see Eq. ( 1271) ]. the absolute values of the effective mass 
and the coupling constant are crucial to the possible (iFSD instability over other instabilities. 
In this sense, microscopic derivation of m and g as well as magnetic field dependences of 
other instabilities are important future issues. 

Some order competing with the rfFSD instability is in fact expected in Sr3Ru207. While 
the experimental phase diagram 17| is very similar to our result Fig. [T](a), the closer com- 
parison between them reveals a difference of slope of the first-order-transition line. In the 
experiment, the edges of the first order line are shifted to the center of the phase diagram so 
that the (iFSD state is stabilized in a narrower region for lower T. This can be interpreted 
as development of some ordering tendency for lower T in Sr3Ru207, which then suppresses 



17, a theoretical 



the (iFSD instability. Although a different interpretation was given in Ref. 
result consistent with this interpretation was indeed obtained in the case of competition of 
the (iFSD and superconductivity. |43l] 

It should be kept in mind that even if the (iFSD instability does not become a leading 
instability, the system can still keep appreciable correlations of the (iFSD.[l^ As a result, 
small external perturbations such as an anisotropic strain and a lattice anisotropy can drive 
sizable (i-wave type FS deformations. This idea was proposed for high-temperature cuprate 



superconductors. [l| 
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Sr3Ru207 is frequently refereed to as a material with a metamagnetic QCEP. 3l| The 
compelling evidence for this was the systematic decrease of a metamagnetic critical end 
point by rotating the magnetic field out of the plane. 
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In fact, several theoretical scenarios 



based on a metamagnetic QCEP were proposed. |35l. I37l| However, recent data for ultrapure 
crystals showed that the critical end point does not reach zero. 40|] Quite recently we have 
found that several important features are already well captured even without the putative 



QCEP. 



heat 
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341 The non- Fermi liquid like behavior of the magnetic susceptibility 201] and specific 



27 



28| observed in Sr3Ru207 can originate from the van Hove singularity of the 



density of states. Moreover, while the first c/FSD transition as a function of a rnagnetic 
field is generically accompanied by a metamagnetic transition [Fig. [T](a) and Ref. [33|, the 
(iFSD transition does not lead to a metamagnetic QCEP in the present model. It remains 
to be studied whether a concept of a metamagnetic QCEP can be a good basis to discuss 
electronic properties in Sr3Ru207 and how the putative QCEP can be related to the (iFSD 
instability and the van Hove singularity. In this sense, it is important to clarify whether the 
anomalous T dependence of the resistivity observed around the metamagnetic transition 261] 
can be explained in terms of c?FSD fluctuations and the van Hove singularity or whether we 
have to invoke quantum fluctuations originating from some QCEP. 

Different scenarios from the rfFSD and the QCEP were proposed for Sr3Ru207, micro- 



scopic phase separation due to Coulomb energy |38l] and magnetic domain formation due to 
long-range dipolar interactions. [39^ In our analysis, we employ the ground canonical ensem- 
ble and a first order phase transition occurs as a function of a magnetic field or the chemical 
potential for low T. Such a transition in turn appears as a magnetic phase separation or 
an electronic phase separation as a function of its canonical conjugate variable, namely the 
magnetization [see Fig. [11(h)] or the electron density [see Fig. Wis) i^^ R-^f. 



10| . Therefore 



invoking the dipolar interaction or Coulomb force, we expect some inhomogeneous states in 



line with Refs. 



38 



and 



39|. However our possible inhomogeneous state may replace only phase 



separated regions in Fig. HTh) and thus is realized only near a metamagnetic transition, in 



contrast to Refs. 



38 



and 



39 



, where the inhomogeneous state is stabilized in the entire region 



of the phase diagram. 
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VI. CONCLUSION 



We have performed a comprehensive analysis of magnetic field effects on the dFSD in- 
stability in a one-band mean-field model with a pure forward scattering interaction. In the 
plane of /i and h{> 0), the (iFSD instability occurs around the axis of /i = ±h, namely 
around the van Hove filling of the majority band (/i < 0) and the minority band (/i > 0), 
even in the limit of h ^ oo. The magnetic field, however, suppresses substantially the energy 
scale of the (iFSD instability. The (iPSD instability occurs through a second order transition 
at high T and typically changes to a first order transition at low T; the end points of the 
second order line are tricritical points. Neither a QCP nor a QCEP of the c/FSD is realized 
by the magnetic field in the present model. In the weak coupling limit, typical quantities 
characterizing the phase diagram have the field- independent single energy scale, while its 
dimensionless coefficient varies with the magnetic field. The magnetic field-induced (iFSD 
instability is a promising scenario for Sr3Ru207 and we have discussed various future issues 
to establish such a scenario. 
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APPENDIX A: PHASE DIAGRAM FOR fi = -0.35 AND -0.34 

Figures [TJ^a) and [2] are two typical phase diagrams as a function of h in the present model. 
These phase diagrams, however, are not connected smoothly by changing fi. In fact, two 
other types of the phase diagram appear in a very limited /i region as shown in Fig. [HI For 
fi = —0.35 [Fig. [8](a)], the dome-shaped transition line is realized on the relatively large 
h side. In addition, the region surrounded by a first order transition line appear around 
{h,T) = (0,0). At T = 0, therefore, there is reentrant behavior as a function of h and the 
symmetric phase appears in the intermediate h region (0.12 ^ h < 0.14). For fi = —0.34 
[Fig. El^b)], the overall shape of the phase diagram is the same as that shown in Fig. [21 but 
a first order line appears for high T near h = 0, accompanied by a tricritical point. As a 
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result, a first order transition happens as a function of T in a sizable h region. This is a very 
special case in our model since a first order transition as a function of T is usually realized 
in a very limited h region as seen in Figs. [D^a) and [2J 

These peculiar types of the phase diagrams are understood from Fig. E^a) or similarly from 
Fig.m^a). The first order transition line in Fig.O^a) is almost straight near ^ LSe^e^^^^'^^^ 
[fi ^ —0.34 in Fig. 111(b)]. Therefore we can have an extended h region of the first order 
transition as seen in Fig. [HJ^b). When we look at closely the region near ^ l.Se^e'^^^'^^^ 
and h < 0.54eAe~^'^^^^^ (inset of Fig. [5]), the first order transition line turns out to have a 
small inward curvature. This is why a symmetric phase is intervened between the two dFSD 
phases in Fig. M^a). 



APPENDIX B: SPIN-DEPENDENT d-WAVE FERMI SURFACE DEFORMA- 
TION 



We have analyzed the d-wave Fermi surface symmetry breaking in the charge channel. 



From the point of view of Landau Fermi liquids, we can consider Fermi s ur: 
also in the spin channel. This possibility was pursued in several references 



45 



ace instability 



46 



47 



the context of a general Landau Fermi liquid theory, 
hidden order in URu2Si2 



46 



47 



48| in 



48l | and a possible relation to a 



45l | Here we clarify the relation between the charge-channel ciFSD 



and the spin-channel (iFSD, and discuss its relevance to Sr3Ru207. 
The minimal model for the spin-dependent (iFSD reads 

^ = E(^k-/^)^k. + ^ J2 /kk'Sk-Sk'-ZiJ^ank., (Bl) 

where Sk = | J2a p'^ka^aisCkis, /^k' = ~5'"c^k'^k'5 and the rest of notation is the same as the 
model ([1]). Since the interaction has SU(2) symmetry and the magnetic field is assumed 
to be applied along the z direction, we assume that the Sz component can have a finite 
expectation value. Defining a mean field 



we obtain the mean-field Hamiltonian 



„a 



N 



(B2) 



El V 



a\2 



k,(T 



(B3) 
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where ^^o- = + 1'^'7^'^k — /^o-- The grand canonical potential per lattice site thus reads 

u; = --y log(l + e-«-/^) + ^ (B4) 

= --y log(l + e-«-/^) + ^ . (B5) 

k,cr ^ 

In the second line, we have introduced 

= 4 + Iv^d^ - f^. , (B6) 

noting that changes its sign with respect to kx ^ ky so that crrfk in the original can 
be written as d^, in Eq. (]B6|) . 

Comparing Eqs. flB5l) and flB6l) with Eqs. ([7]) and ([6]), respectively, we see that the free 
energy becomes exactly the same under the following mapping, 

Iv^^V, \9'^9- (B7) 

Hence the thermodynamics in the spin channel of the dPSD is the same as that in the charge 
channel in the sense that we obtain the exactly the same results as Figs. W^sl) and (e)-(h) 
under the mapping (1B7I) . The difference appears in the "internal" structure of the order 
parameter and in a deformation of the FS. In the spin channel, we can write |?7" = | cr?7^, 
where 



(B9) 



2 AN 

k 



-SEk^k/leg't) fora=T 
^|jEk4/(e£l) fora=i . 



Comparing with Eq. (Q, we see the relation under the mapping ( \B7\) 

Iv'^^Vu —Vl^Vi- (BIO) 

As we have actually seen in the numerical result in Sec. Ill, ?7| in general has the same sign 
as 7]^ in the charge channel. Therefore 1]'^ has an opposite sign of r/" as seen in Fig. M^a). 
While we have introduced [Eq. (1B6P ]. the Fermi surface itself is defined by = 0. 
Since contains a factor ar]°'dii_, a Fermi surface deformation in the spin channel occurs 
always in the opposite direction between the up-spin and the down-spin electron band as 
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shown in Fig. Mjo); note that the deformation is determined by 77", not by t]^. As a result, 
the net deformation of the band is partially compensated. This is a crucial difference from 
the (iFSD instability in the charge channel. 



The recent experiment by Borzi et al. 



32| showed a strong a;?/-anisotropy of the magne- 



toresistivity by applying an additional small magnetic field to the RUO2 plane. This strong 
anisotropy may be discussed more naturally in terms of the (iFSD instability in the charge 
channel rather than the spin channel. A conclusive discussion on which channel is more 
dominant would be to study microscopic deviation of the c/FSD attractive interaction in 
both charge and spin channel in the context of Sr3Ru207 and to compare the strength of 
the each channel. 
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FIG. 1: The mean-field solution for t' /t = 0.35, 5 = 1, and /x = —0.4. (a) h-T phase diagram; 
the transition line contains a second order line, T^"*^, for high T and two first order line, T^^^*, for 
low T; the solid circles are tricritical points; the dotted line denotes the van Hove energy of the 
up-spin band, (b) h dependence of the order parameter at T = 0.01; note that r] = -q^ + r^. (c) 
and (d) FSs for h = 0.3 and 0.5 at T = 0.01; the solid lines (gray lines) are FSs for g = 1 ((7 = 0); 
the deformation of the inner FS in (d) is hardly visible, (e) h dependence of n^- at T = 0.01 for 
g = 1 (solid line) and (dotted line), (f) h dependence of the magnetization at T = 0.01 for ^ = 1 
(solid line) and (dotted line), (g) T dependence of the magnetization for several choices of h\ the 
dotted line represents the result for 5 = 0. (h) m-T phase diagram; the system undergoes phase 
separation in the shaded regions surrounded by T^^; the other notation is the same as that in (a). 
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FIG. 2: h-T phase diagram for t' /t = 0.35, 5 = 1, and ^ = —0.2; the notation is the same as that 
in Fig. d^a). 
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FIG. 3: h dependence of m at ^ = —0.2 and T = 0.01 for 5 = 1 (sohd line) and (dotted hne). 
The inset shows the h dependence of the order parameter. 
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FIG. 4: (Color online) (a) The first order dFSD transition line in the plane of /u and h at T = 0.01 
for g = 1; the symmetry-broken phase is stabilized in the colored (shaded) area; for notation, see 
the text, (b) and (c) h dependence of Tc along the van Hove energy, namely along the dotted line 
in (a) for /i < and /i > 0, respectively; the solid (dotted) line is the result for = 1 (0.5); T^^ 
for g = 0.5 is multiplied by 6. 




FIG. 5: (Color online) The first order dFSD transition line in the weak coupling limit in the plane 
of /X and h at T = 0; /x and h are scaled by the energy eAe"^''^^^-'; the symmetry-broken phase 
is stabilized in the colored (shaded) area; the dotted line represents the van Hove energy of the 
up-spin band (// < 0) and the down-spin band (fi > 0). The inset magnifies the region around 
fi « 1.3. 
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/z/r/"(0) 

FIG. 6: Tc along the van Hove energy, namely along the line of /i = |/i|, in the weak coupling limit; 
Tc and h are scaled by TJ^{<d). 
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FIG. 7: (Color online) (a) /x dependence of , T*''\ /itri — h^n, and /ii — /ivH, where /ivH = 
1^1; all quantities are scaled by the energy e\e~^/^'^^\ (b) /i dependence of the universal ratios 
T*"/l^tri-/ivH| andrt"/r,-H. 
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FIG. 8: h-T phase diagram for fi = —0.35 (a) and fi = —0.34 (b) for t'/t = 0.35 and g = I; the 
notation is the same as that in Fig. HJa). 




FIG. 9: The mean-field solution in the spin-channel dFSD instability for t'/t = 0.35, fx = —0.4, 
and g"" = 4. (a) h dependence of the order parameter at T = 0.01. (b) FSs at /i = 0.3 (solid lines); 
the gray lines denote the FSs for g^ = 0. 
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